Optimal generation of entanglement under local control 
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We study the optimal generation of entanglement between two qubits subject to local unitary 
control. With the only assumptions of linear control and unitary dynamics, by means of a numerical 
protocol based on the variational approach (Pontryagin's Minimum Principle), we evaluate the 
optimal control strategy leading to the maximal achievable entanglement in an arbitrary interaction 
time, taking into account the energy cost associated to the controls. In our model we can arbitrarily 
choose the relative weight between a large entanglement and a small energy cost. 
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Introduction. — The building block for the implementa- 
tion of quantum technologies is a pair of interacting two- 
level systems (qubits), whose evolution can be affected 
by external actions. The basic operations to be accom- 
plished through this system are the storage and manipu- 
lation of encoded information. The use of this prototype 
of physical apparatus is motivated by its quantum nature, 
that is the existence of some peculiar properties not ex- 
hibited by classical systems. Among them, the quantum 
correlation called entanglement plays a prominent role in 
the realization of outperforming protocols. 

The standard setting for the manipulation of a pair of 
qubits is the so-called Local Unitary control [H, HJ. In 
this framework, neglecting the influence of the external 
environment, the state system is represented by the unit 
norm complex vector \ip(t)), satisfying the Schrodinger 
equation 



\m) = -iH T ( ul (t),u 2 (t))\ip{t)), 



(i) 



where H T ( Ul {t), u 2 (t)) = + H 2 (u 2 (t)) + Hi is 

the total Hamiltonian. The local contributions Hi(ui(t)) 
(i = 1,2) can be modified by means of external actions, 
represented by the control functions U\ (t) and u 2 (t) , and 
Hi is the uncontrollable interaction term, responsible for 
the entanglement created in the system. We find conve- 
nient to express the initial state | ^0) = |"0(O)) according 
to the Schmidt decomposition, 



I Vo) = VP <8 |x) + y/i-P \<p) x ® lx>" 



(2) 



where P € [0, 1] and (<p\ip) = (x|x) = 0- K there is not 
initial correlation between the two qubits, P = or P = 1 
and jV'o) is a product state. Conversely, a maximally 
entangled state (Bell state) corresponds to P = 0.5. As 
a measure of entanglement we introduce the concurrence 
C denned as 



c\t) = \m)Wy®v y m))*\ 



(3) 



assuming values in the interval [0, 1], vanishing for uncor- 
rected states, and reaching its maximum for maximally 



entangled states. This quantity satisfies all the properties 
of an entanglement monotone Q . 

Because of it fundamental relevance, several issues re- 
garding the entanglement generation in the system ([I} 
have been addressed in the past years, as well as the 
problem of generation of nonlocal gates. In Q, Dur et 
al. characterized the capability of creating entanglement 
for a generic interaction Hi. In particular, they provided 
a strategy to minimize the time of generation of entan- 
glement through arbitrarily fast local control, and con- 
sidered the impact of ancillas. Moreover, these authors 
proved that an initial amount of entanglement leads to a 
more efficient production of entanglement. In a different 
context, Kraus and Cirac expressed the maximal attain- 
able entanglement for an arbitrary unitary operator, as 
well as the corresponding initial factorized state [Bj . The 
generation of entanglement implemented via local mea- 
surements has been considered in [f| . 

Timc-optimality has been further considered in the 
context of the simulation of a quantum gate, by intro- 
ducing the Interaction Cost, that is the minimal time to 
perform a gate using local operations 0, B 0] • 

While time-optimal procedures are fundamental for 
the implementation of efficient computational architec- 
tures, they usually ask for impulsive controls (instan- 
taneous local manipulations of arbitrary strength), and 
their operational cost (the energy loss associated to the 
local controls) is not accounted for. Moreover, for some 
physical apparatus it could be difficult to precisely set 
the optimal interaction time, and a predetermined time 
could be preferable. 

With these motivations in mind, in this letter we de- 
scribe a control theoretical approach accounting for the 
aforementioned cost and considering an arbitrary inter- 
action time. Under the assumption of linear control with 
Hi(Ui(t)j = Ui(t) Hi, i — 1,2, and using a numerical 
protocol based on the variational method, we evaluate 
the best control strategy, that is the optimal control 
functions Ui(t) driving an arbitrary initial state |^>o) &s 
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close as possible to a maximally entangled state. This is 
not the most general form for linear control, however it 
highly reduces the computational complexity of the prob- 
lem while preserving most of its relevant features. This 
approach represents a novelty with respect to previous 
treatments and it complements them. We will mainly 
refer to 0] for a comparison. 

Optimal control methods have been initially used in 
quantum mechanics for the control of molecular dynam- 
ics [Toj j . Recent applications of these techniques in quan- 
tum information are in particular aimed to determine 
the optimal gate generation [ll|, 12, 13|and the optimal 
evolution and state transfer (e.g. see 3, la])- 

The computational procedure. — We find the optimal 
control strategies Ui(t) using an iterative procedure based 
on the variational approach known as Pontryagin Mini- 
mum Principle. The computational tools employed to de- 
rive the results presented in this work can be used as well 
for the solution of more general optimal control problems. 
In fact, different performance measures as well as more 
complicate system dynamics (for example in the presence 
of irreversibility and dissipation) can be imposed. A de- 
tailed discussion of the protocol, of its performance, and 
of further applications is out of the scope of this letter, 
and it will be presented in a forthcoming paper. How- 
ever, for sake of completeness, we summarize here the 
basic ideas underlying the procedure, without entering 
into details. 

Consider a system described by the state x(i) : E — > 
E™ whose dynamics and initial conditions are given by 



x(t)=f(x(i),u(t)), x(0)=x o , 



(4) 



where u(<) : E — ► M m is a vector of control functions and 
f(x, u) : E n x W m -> E n is a vector field. Assume that 
u(t) has to be chosen such that the cost functional (or 
performance measure) 

J(u(t)) =$(x(r)) +a [ jSf(x(t),u(f)) dt (5) 
Jo 

is minimal, where r is the fixed final time and $(x) : 
W 1 -> E, (x, u) : E™ x R m -> E are arbitrary functions. 
Notice that the cost J contains a final-time term and an 
integral contribution (that we shall later denote by I(t)), 
and the real coefficient a expresses the relative weight of 
them. Define the optimal control Hamiltonian as 

Jf( X (t), p(t), u(t)) = f(x(t), u(i)) • p(t) + JSf (x(t), u(t)) , 

(6) 

where • is the inner product, and p(t) : E — > E™ is an 
auxiliary variable conjugate to x(t), often called costate, 
whose dynamics is 



p(t) = -V x ^(x(t),p(t),u(t)) 
with final condition 

P(r) = -V x $(x(r)). 



(7) 



(8) 



Then the Pontryagin Minimum Principle can be stated 
as follows. 

Proposition If we denote by u'(t) the optimal control 
strategy, and by x'(i) and p'(i) the corresponding opti- 
mal state and costate trajectories, then, for all t £ [0, r], 

Jf( x '(t),p'(t),u'(t)) < ^(x'(t),p'(t),u(t)). (9) 

Consequently, 

v u jr(x / (t), P / (t),u'(t)) =o. (io) 

While a more general formulation can be given to this 
principle (in particular a non-fixed final time r can be 
considered), this approach is all we need for our purposes. 
For more details and for the proof of the principle, see 
[l6| and [13] (mainly focusing on quantum mechanical 
applications). 

In our case, we find convenient to represent the 
state vector in the computational basis {\ei(t)),i — 
1, . . . , 4} given by tensor products of eigenvectors of <7 2 , 
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1=1 



(11) 



and x(i) = (Reipi(t),hnipi(t);i = 1,...,4) , where T 

means transposition. Moreover, u(t) — (ui(t), U2(t)) T , 
and the vector field f is linear in both x and u. The final- 
time contribution to the cost function is the deviation of 
the final entanglement from its maximal attainable value, 



d>(x(r)) = 1 - C(t) 



(12) 



whereas the integral part measures the energy loss, as- 
sumed to be proportional to the squared norm of u, 



J2?(u(t)) =ul(t)+u 2 2 (t), 



(13) 



the prototype of energy cost for a nuclear spin driven by 
a magnetic field. 

Our protocol consists of an iteration in which, after 
solving (|3]) and with boundary condition ((H)), the con- 
trols are redefined step by step in order to fit the condi- 
tion (|10j) . The procedure starts with arbitrary trial func- 
tions Ui(t) and it stops when a predetermined accuracy 
level is reached. 

Discussion of numerical results. — The optimal time 
strategy described in [J] is based on the maximization 
of the entanglement rate at every time. This procedure 
leads to a vector \i/j(t)) whose (time-dependent) Schmidt 
coefficient is given by 



P(t) =sin 2 (h max t + (j) Q ), 



(14) 



where P(0) = P = sin 2 0o- The entanglement capability 
hmax measures the ability of the interaction to produce 
entanglement, its definition and expression in terms of 
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FIG. 1: Dependence on r of the final entanglement C(t) and the 
energy cost I(t) for several values of the weight parameter a. The 
initial state is an optimal state for the Heisenberg Hamiltonian, 
with P = 0,\<p) = \ T)z, and |x) = | 



FIG. 2: Dependence on r of the final entanglement C(t) and the 
energy cost J(t) for several values of the merit function a. The 
initial state is given by P = 0, |^>) = | [) y , and |x) = | i)z- 



the singular values of Hi are given in [4J. The opti- 
mal time r op t is defined as the smallest time such that 
P(t) = 0.5. Equation (|14| defines the steepest entangle- 
ment growth, therefore it represents the upper limit for 
the production of entanglement up to the optimal time. 

Our protocol works for arbitrary Hamiltonian terms 
and initial states. In order to illustrate the main results 
of this work we need to fix them. Because of its relevance, 
we consider the Heisenberg interaction Hj — a x ® a x + 
o y (g> <jy + a z <g> a z (with h max = 2 and T opt = n/8), 
Hi = <j x ® I, and H2 = I (E> <J X . 

In Fig. [T] and [2] the dependence on r of the final en- 
tanglement C(t) and the energy cost I(t) (the integral 
in |5])) is shown for several values of a, for the optimal 
strategies. The two figures correspond to different un- 
corrected initial states (P — 0): \<p) = \ f) z , |x) = | \) z 
and \<p) = I i) y , \x) = I i) z respectively (| |)i and | j), 
denote the +1 and —1 eigenvectors of the Pauli matrix 
en, i = x,y,z). 

The first state is an optimal initial state for the Heisen- 
berg interaction 0, @, then it evolves in the optimal 
time to a maximally entangled state without local ac- 
tions. This is apparent from Fig. [TJ all the represented 
curves for C(r) have the same behavior for r ^ r opt , 
independently of a. These patterns fit the entanglement 
evolution associated to (fT4]) . represented by the grey line, 
since the system is driven by Hi along the optimal tra- 
jectory. Therefore, the energy cost associated to these 
paths vanishes. For r > T opt the parameter a becomes 
relevant, since a local action is necessary to maximize the 
final entanglement. We observe that the goal C(r) = 1 is 



approached as a is decreased, that is an higher cost is tol- 
erated, and that, for every a, C(r) — > 1 when r increases. 
The peaks of I(t) correspond to the valleys associated to 
(fT4| . where a stronger local control is needed, and the 
magnitude of these peaks decreases with r. 

In Fig. [2] a similar analysis is presented for a particu- 
lar non-optimal state. In this case, local manipulations 
are needed even for r ^ r opt since the state has to be 
adapted in order to fully exploit the entangling capabil- 
ity of the interaction. This is apparent in the presented 
plots, where the growth of C(t) is steeper with a smaller, 
and eventually it approaches the optimal curve associated 
to (|14| Correspondingly, there is a relevant initial contri- 
bution to I(t). Notice that for an arbitrary non-optimal 
state, different patterns could be found for a large, in 

u (t) 




FIG. 3: Optimal strategies u\(t) and U2{i), and corresponding 
entanglement C(t), for the optimal initial state, and r = 4. We 
have chosen a = 10 — 1 . 
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FIG. 4: Dependence on r of the final entanglement C(r) and the 
energy cost I(t) for several values of the initial Schmidt coefficient 
P. The initial state is defined by \tp) = | f) z , and |x) = I l)z- We 
have chosen a = 10~ 2 . 



particular the steepest growth of C(r) could be slower 
than the optimal one, and the plateau of C(r) could be 
(even significantly) below 1 . This is due to the particular 
choice of Hi and Hi , that could be inappropriate for the 
initial state considered. 

The optimal control strategy can be fixed in line with 
these considerations. Given the interaction time r, the 
choice of the coefficient a represents a compromise be- 
tween magnification of C(r) and reduction of I(r). The 
optimal strategies u\ (t) and ui (t) match the requests on 
r, C(t) and I(t). An example is provided in Fig. [3J 
with the Hamiltonian terms and the optimal initial state 
previously introduced, and r = 4. 

Oscillating controls with modulated amplitude are usu- 
ally obtained when r exceeds the optimal time. When the 
energy cost is a relevant factor, impulsive controls are 
less efficient than controls distributed over time. They 
become the optimal strategy when r < r opt , if a large 
energy expense is accepted. The function C(t) does not 
in general stabilize around its maximum (as in the exam- 
ple presented here) , however the amplitude of oscillations 
usually decreases as t — > r. 

We have also considered initial states with a non- 
vanishing entanglement. Some plots are shown in Fig. 
Hwith \<p) = | t) 2 , |x) = | i)z and P ^ 0. As intuition 
suggests, the cost I(r) is usually smaller for correlated 
initial states. 

Conclusions. — We have considered a variational ap- 
proach for the solution of optimal control problems in- 
volving two qubits. By accounting for the energy cost 
associated to the manipulations of the system, we are 
able to find the optimal strategies to be used in order to 
drive the system. 

In this letter, we have described the entanglement gen- 
eration for systems without interaction with the external 
environment, driven by local unitary control. Our numer- 
ical analysis is consistent with previous results, in par- 
ticular it reproduces the optimal entanglement growth 
and the corresponding minimal time. Moreover, in our 
approach the interaction time can be arbitrarily chosen 



(for example, it can be a fixed instrumental time or a 
predetermined operational time), it is not fixed by the 
interaction. From this point of view, our protocol com- 
plements the existing methods for the generation of en- 
tanglement, and it is of interest whenever a non-optimal 
interaction time is preferred or the energy cost associated 
to the controls has to be taken into account. 

Using standard Hamiltonian terms and particular ini- 
tial states, we have studied the relations among the rel- 
evant quantities, and provided some examples. We have 
found that a large interaction time is usually preferred 
for the reduction of the energy cost, without decreasing 
the efficiency of the entanglement production. From this 
point of view, controls spread over time are more conve- 
nient with respect to impulsive controls. 

R. Romano acknowledges support from the European 
grant ERG:044941-STOCH-EQ. Work in part supported 
by INFN, Sezione di Trieste, Italy. 
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